The microcanonical temperature in the Wang-Landau sampling 
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In this paper we present a brief comment on the behavior of the slope of the entropy in the ground 
state. We show that the inverse temperature of the ground state diverges only for an infinite system 
and therefore the acute angle that the slope displays at this temperature is related to finite size 
effects. 



Wang-Landau sampling ( l/KLS')[ll, ^ is the most 
widespread technique among the entropic algorithms. It 
simulates efficiently the density of states in a given test 
model. The heuristic idea of the method is based on the 
fact that if one performs a random walk in energy space 
with a probability proportional to the reciprocal of the 
density of states, a flat histogram is generated for the 
energy distribution. Since the density of states produces 
huge numbers, instead of estimating g{E), the simulation 
is performed for S{E) = Ing(iJ), and a histogram H{E) 
is accumulated during the simulations to control the fre- 
quency of visits to the energy levels. At the beginning of 
the simulation we set S{E) = for all energy levels. The 
random walk in the energy space runs through all energy 
levels from Emin to E^ax with a probability 

p{E E') = min(cxp [S{E) - S{E')],1), (1) 

where E and E' are the energies of the current and the 
new possible configurations. For each given model one 
should define a Monte Carlo sweep (e.g. spin-flip trial 
for the 2D Ising model or N monomer moves for a ho- 
mopolymcr). After each Monte Carlo sweep we update 
H{E') H{E') + 1 and S{E') S{E') + F, where F = 
In / (/ is the so-called modification factor) . The initial 
modification factor is taken as / = /o = e = 2.71828.... 
The flatness of the histogram is checked after a number 
of Monte Carlo (MC) steps and usually the histogram is 
considered flat if H{E) > 0.8(H), for all energies, where 
(H) is an average over the energies. If the flatness con- 
dition is fulfilled we update the modification factor to a 
finer one by setting /i+i = ^/Ji and reset the histogram 
H{E) ~ 0. Testing the behavior of the canonical aver- 
ages during the simulations, one can determine the ideal 
modification factor to stop the simulations Q . Having in 
hand the density of states, one can calculate the canoni- 
cal average of any thermodynamic variable as 



(2) 



where (X) e is the microcanonical average accumulated 
during the simulations and /3 = l/Zc^T, kg is the Boltz- 
mann constant and T is the temperature. One of the 
interesting features of the method is that it can also ac- 
cess some quantities, such as the free energy and entropy 
at any temperature, which are not directly available from 
conventional Monte Carlo simulations. 



However, analyzing the simulation results in terms of 
the microcanonical ensemble, the logarithm of the den- 
sity of states corresponds exactly to the definition of en- 
tropy 



S^kg hig{E). 



(3) 



The Boltzmann constant ks ensures the agreement with 
the Kelvin scale of temperature, defined by 

j_ _ as; 

T ~ dE' 



(4) 




FIG. 1. Density of states obtained by Wang-Landau sampling 
for the 2D Ising model for L = 32, 64. 

In this work we intend to shed some light on the reason 
why the slope of the density of states is not infinite in the 
ground state, where the inverse temperature diverges. In 
Fig. [1] we show the logarithm of the density of states of 
the 2D Ising model obtained by Wang-Landau sampling 
for two lattice sizes: L = 32, 64. One can see that indeed 
the slopes are not infinite in the ground state. More- 
over, this behavior is also observed in the exact results 
of BealeQ. The reason for this apparent contradiction 
are the finite-size effects. This situation resembles that 
of the magnetization in the Ising model, which for an 
infinite system drops to zero exactly at the critical tem- 
perature but presents a tail for finite sizes. 

For any finite lattice in the 2D Ising model, the smallest 
AE from the ground state is given by 



AE = {Er, 



En 
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FIG. 2. Finite-size scaling behavior of -44 as function of L. 



Therefore, the limit in the derivative 

AS 



dS_ 

dE 



hm , „ 

AE^O AE 
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becomes exact only if L — > oo {Emin = — — >■ — oo). 

In order to investigate the finite-size scaling be- 
havior of the ratio AS/AE we performed Wang- 
Landau simulations of the 2D Ising model with L = 
32,36,40,44,48,52,56,64,72 and 80, taking N = 
24,24,20,20,20,16,16,16,12 and 12 independent runs 
for each size, respectively. 

Setting ks equal to unit, S{E) = \og{g{E)). In FigH 
we show that the derivative in the ground state scales as 



dS{L) 
dE 



= AL", 



(7) 



where A = 0.3605(44J and a = 0.2546(31). 

One can see that does actually tends to infinity as 
L — ^ oo. The misunderstanding is therefore related to 
finite size effects. 
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